Inhomogeneous level splitting in Pr 2 - x Bi z Ru 2 07 
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We report that Bi doping drives Pr2-nBi a; Ru207 from an antiferromagnetic insulator (a; = 0) to 
a metallic paramagnet (x ~ 1) with a broad low T maximum in C/T. Neutron scattering reveals 
local low energy spin excitations (Huj ~ 1 meV) with a spectrum that is unaffected by heating to 
fcsT >> hut. We show that a continuous distribution of splittings of the non-Kramers Pr 3+ ground 
state doublet such as might result from various types of lattice strain can account for all the data. 



PACS numbers: 75.10.Dg, 71.70.Fk, 71.27.+a 

Rare earth ions with an even number of electrons 
in a partially filled 4f-shell require a highly symmetric 
environment to retain ground state degeneracy. Such 
non-Kramers degeneracy is consequently unstable to- 
wards lattice distortions, which can occur spontaneously 
through the Jahn- Teller effect, dynamically through in- 
teractions with phonons, Q or at random due to static 
inhomogeneities 0. In this letter we describe an in- 
advertent encounter with the latter effect while explor- 
ing the unusual thermo-magnetic properties of a diluted 
frustrated magnet, P^-^Bi^R^Oy. While the enhanced 
specific heat of the material is reminiscent of ultra heavy 
fermion behavior, the temperature (T) and wave vector 
(Q) independent neutron spectra and the diverging sus- 
ceptibility implicate inhomogeneous level splitting of the 
non-Kramers Pr doublet. Our comprehensive analysis 
links neutron, specific heat, and susceptibility data and 
is important because similar physics may be operative in 
other magnets of recent notoriety. 

The parent material, Pr2Ru20y, contains two sub- 
lattices of magnetic ions on the vertices of corner-sharing 
tetrahedra. This lattice is notable because antiferromag- 
netic nearest neighbor interactions do not select a unique 
long range ordered spin configuration but a highly degen- 
erate manifold. While Pr 2 Ru207 does develop magnetic 
order at T/v = 165 K, dilution of the Pr 3+ sites by the 
smaller isovalent Bi 3+ ion converts the ordered insulator 
to a metallic paramagnet. Similar behavior was found in 
other rare earth ruthenates, however, the Pr based ma- 
terial has a greatly enhanced low— T specific heat with 
analogies to non-fermi- liquid and heavy fermion systems. 
Here we link the anomalies to static inhomogeneous split- 
ting of the non-Kramers Pr 3+ doublet, and consider var- 
ious options to account for the inhomogeneity. 

Powdered samples of P^-^BicRi^Or were synthesized 
by the solid state reaction method. The mixtures of 
P1-2O3, Bi2C>3 and RuC>2 in proper molar ratios were 
pre-reacted at 850 °C for 15 h in air and then ground 



and pressed into pellets. The pellets were sintered at 
1000-1200 °C in air with intermediate grindings. Powder 
X-ray diffraction measurements indicated single phase 
samples. Powder neutron diffraction data were collected 
on a x ~ 1.0 sample, using the BT1 diffractometer at 
NIST. Rictveld analysis confirmed the cubic pyrochlorc 
structure and inductively coupled plasma-emission spec- 
troscopy showed the Pr:Bi composition to be 1:0.94 cor- 
responding to x = 0.97. This sample was subsequently 
used for inelastic neutron scattering experiments. 

Zero field heat capacity measurements were carried out 
on all members of the series using a commercial calorime- 
ter. Field dependent heat capacity measurements cm- 
ployed a home made probe with a Si platform in a su- 
perconducting magnet. The electrical resistivity, p, of 
the polycrystallinc samples was measured by the four- 
probe method. The overall temperature dependence of 
p for different doping levels, x, are similar to those re- 
ported for Y 2 _ a: Bi a; Ru207 0,11 Pr 2 Ru 2 7 {x = 0) is 
insulating while Bi2Ru2C>7 [x = 2) is metallic (Fig. QJ. 
For intermediate values of x, p decreases abruptly with 
x indicating an insulator to metal transition for x « 0.8. 

Magnetic susceptibility data were collected on selected 
samples using a SQUID magnetometer. The insulating 
x = sample shows abrupt development of hysteresis for 
T < Tn = 165 K indicating Neel order plus spin canting 
from anisotropic exchange interactions pi . While the 
x = 2 sample behaves as a Pauli paramagnet Q , samples 
near the metal insulator transition for x ~ 0.8 have an 
increasing susceptibility upon cooling but fail to develop 
magnetic order down to T = 1.5 K (Fig.[5J). 

Fig-Ea) shows that metallic Pr2- x Bia;Ru2 07 samples 
close to the metal insulator transition have a strongly en- 
hanced low— T specific heat, C(T). In this figure the lat- 
tice contributions have been subtracted using data from 
the x = 2 sample corrected for unit cell mass through the 
Debye interpolation scheme. For x = 0.8 and T « 2.0 K 
the specific heat (C/T = 0.7 J/mole-f.u./K 2 ) is more 
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FIG. 2: Temperature dependence of the magnetic suscepti- 
bility x °f Pr2-a;Bi :E Ru2 07 for various x. The solid lines 
show the calculated magnetic susceptibility defined as x(T) = 
X'(T) + xcf(T). Here X '{T) is given by Eq. © and Xcf(T) 
is the contribution from the higher energy CF levels. The 
dashed line shows the calculated susceptibility for a doublet 
ground state. Inset, the Pr 3+ CF energy level scheme. 



FIG. 1: Temperature dependence of the specific heats C/T 
of Pr2- a: Bi a ;Ru207 (a) for various x and H = and (b) for 
various fields (H) for x — 1.0. The lattice contribution was 
subtracted as described in the text. The solid lines show 
C/T calculated from the observed neutron spectra and the 
assumption of a quenched level distribution (Eq. (3J). A 7 
term of 0.04 J/mole/K 2 was added to the model in (b). Inset, 
electrical resestivity, p, of various x at 30 K. 



than three orders of magnitude greater than for a con- 
ventional metal. Figure IHb) shows that the enhanced 
specific heat is very sensitive to a magnetic field indicat- 
ing a high density of magnetic states. 

To determine their origin we performed inelastic neu- 
tron scattering experiments on the HET and IRIS spec- 
trometers at the ISIS Facility. On HET we used inci- 
dent energies of Ej= 35 meV and 160 meV, which pro- 
duced an elastic energy resolution of 1.4 meV and 7 
meV respectively. For the IRIS experiment disk chop- 
pers selected an incident energy spectrum from 1.3 meV 
to 4.6 meV pulsed at 25 Hz and a backscattering py- 
rolytic graphite analyser system with a 25 K Be filter 
selected a final energy of Ef= 1.847 meV for an elas- 
tic energy resolution of 17.5 fieV. Count rates were nor- 
malised to the incoherent scattering from a vanadium 
slab for the HET experiment, while IRIS data were nor- 
malized using coherent nuclear scattering from the sam- 
ple. These procedures yield absolute measurements of 
I{Q,Tuj) = (ir Q ) 2 \£F(Q)\ 2 2S{Q,huj) to an overall scale 
accuracy of 20 %. Here 7= -1.913 and g are spectro- 
scopic g-factors of the neutron and magnetic ion respec- 
tively, ro= 2.82 fm is the classical electron radius, F(Q) 
is the magnetic form factor for Pr Q and S(Q, Tuo) is the 
spherically averaged dynamic correlation function 0, Q . 



To determine the relevant low energy degrees of free- 
dom we used inelastic neutron scattering to map the crys- 
tal field (CF) level scheme for Pr 3+ in the sample. Five 
magnetic excitations were identified at 9.5(1) meV, 50(5) 
meV, 84.0(1) meV, 103.1(1) meV, and 115.7(2) meV 0. 
The number of excitations is consistent with a J = 4 ion 
in D 3 d symmetry. The fact that excitations are observed 
to all 5 excited states implies a doublet ground-state. 
Further analysis of the spectrum indicates a ground state 
wave function of the form |±) = a\ =f4) +j3\ ± 2} ±7! ^fl), 
where a = 0.90(4), (3 = -0.10(3), and 7 = 0.42(7). CF 
excitations in the bismuth doped sample were broadened 
considerably as compared to pure Pr2Ru2 07 where the 
excitations are almost resolution limited. The average 
doping induced half width at half maximum of 3.8 meV 
is a first indication of a distribution of CF environments. 

The low energy magnetic response in the x = 1.0 sam- 
ple features a broad and apparently dispersion-less peak 
centered at ~1 meV energy transfer (Figs. [S] and HJ. The 
large spectral weight and the low energy scale suggests 
that these excitations are related to the enhanced low— T 
specific heat of the material. The wave vector depen- 
dence of the intensity follows the squared single ion form 
factor of Pr 3+ (Fig. |3J), indicating an absence of corre- 
lations between Pr sites and single ion physics. This is 
unusual for a dense Kondo (lattice) system, as in general 
competition between inter-site RKKY and on-site Kondo 
singlet formation produces short range antiferromagnetic 
correlations. The data also appear to be inconsistent 
with isolated Kondo singlets, which cannot scatter inclas- 
tically as Q — > Most revealing, however, 

is that the spectrum, when appropriately scaled (Fig.0J, 
is utterly unaffected by heating to temperatures of order 
the Neel temperature for the parent insulator. This un- 
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FIG. 3: Q— dependence of energy integrated inelastic neutron 
scattering intensity for T = 1.5 K, 30 K, 70 K and 90 K. Data 
were obtained by integrating over 0.5 < huj < 1.5 meV. The 
solid line shows \F(Q)\ 2 for Pr 3+ scaled to the data 0- 
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FIG. 4: Energy dependence of Q— averaged neutron scatter- 
ing data scaled for proportionality to the level distribution 
function, p(A) of Eq. The data collapse indicates that 
p(A) is T— independent in the range probed. 



usual result indicates that the finite width of the sample 
averaged spectrum results from some form of quenched 
inhomogcneity rather than relaxation effects. 

We now pursue such an explanation for low energy 
physics in Pr^-xBi^R^Oy more carefully As Pr 3+ is 
a non-Kramers ion a doublet ground-state relies on D^d 
point group symmetry. Lowering of the CF symmetry 
with an amplitude that varies through the sample would 
split the doublet by varying amounts on different sites. 
Such quenched inhomogcneity would produce a T inde- 
pendent spectrum of singlet-singlet transitions. As it is 
based on single site CF dynamics, this explanation links 
neutron scattering, specific heat, and susceptibility data 
in a parameter free fashion that we shall put to the test. 

Denote the quenched distribution of singlct-singlct 
splitting by p(A) such that Np(A)dA is the number of Pr 
ions with level splitting between A and A + dA and N is 
the total number of formula units in the sample. Using 
expressions for the generalized susceptibility of the CF 
split doublet and the fluctuation dissipation theorem || 
it is straightforward to show that the splitting distribu- 
tion function is related to the sample averaged response 
functions as follows 
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Here j3 = l/k B T, \" denotes the imaginary part of 
the generalized susceptibility, g= 0.8 is the spectro- 
scopic g- factor for Pr 3+ , a 2 is the powder averaged 
squared singlet-singlet matrix element [Dj, and Z'{fi) = 
^2 i riie~P Ei is the partition sum excluding the low en- 
ergy doublet. Fig. H shows {g/2) 2 S(huj){l + e~ 0Hu + 
e -ph<*>/2Z'{[3)), which according to Eq. JIJ and the un- 
derlying assumptions should be T independent. The ex- 



cellent data collapse provides strong evidence that the 
observed spectrum indeed results from inhomogeneous 
splitting of the Pr CF doublet. While the discrepancy 
at negative energies results from incomplete background 
subtraction, the differences in the data for Tilj « 0.1 meV 
can be ascribed to incoherent inelastic phonon scattering 
on account of the approximately linear T dependence. 

Using perturbation theory and Eq. Q the T depen- 
dence of the specific heat and magnetic susceptibility 
can now be obtained directly from the neutron results 
for p(A) with no adjustable parameters. 
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Here A(H Z ) 2 = (2gp B H z (+\J z \+)) 2 + A 2 and H z is the 
component of the magnetic field along the local 3-fold 
axis. Note that Eq. Q becomes the exact Kramers- 
Kronig relation when p(A) from Eq. JIJ is evaluated at 
the running temperature. However, the assumption of 
a quenched distribution backed by Fig. 01 allows calcula- 
tion of the full temperature dependence of x'(^) based 
on a single neutron spectrum. To minimize contributions 
from phonon scattering we used T = 1.5 K data to de- 
rive p(A). Note however that indistinguishable results 
are obtained using the average of all temperatures. 

The lines in Fig. 2] were calculated from Eq. J^J nc ~ 
glecting any variation of p(A) with Pr concentration. 
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We added the nuclear specific heat as calculated from 
the natural isotope distribution of the sample [l4|. In 
the bottom frame we added a common best fit Sommer- 
feld term, 7r u =40 mJ/mole/K 2 , to account for contri- 
butions from the metallic state. This value is not dis- 
similar from 7r u =72 mJ/mole/K 2 for Pauli paramag- 
netic YBiRu 2 C>7. The magnetic susceptibility was calcu- 
lated as x(T) = x'(T) + XCf(T), where X '(T) is given 
by Eq. @J and Xcf{T) accounts for contribution from 
higher energy CF levels [J^. The discrepancy between 
the lines and the data points is of order the accuracy 
in the overall normalization of the specific heat, suscep- 
tibility, and neutron data. As there arc no adjustable 
parameters this supports the inference that Pr ions in 
Pr 2 _ x Bi x Ru 2 07 form a distributed two level system. 

The measured distribution function provides im- 
portant clues to the origin of level splitting in 
Pr 2 _ x Bi x Ru 2 C>7. The absence of T dependence to p(A) 
from 1.5 K to 90 K (see Fig.EJl all but rules out dynamic 
phenomena. Fluctuations slow enough to be quasi-static 
on the time scale of the « 1 meV singlet-singlet tran- 
sitions would be thermally activated in that tempera- 
ture range. So rather than reflecting changes with time 
for a single Pr ion, p(A) represents a spatial distribu- 
tion of Pr environments. Because Eq. @ accounts for 
finite field data too, magneto-static disorder from high T 
Ru spin freezing seems unlikely though this deserves fur- 
ther examination through pSK measurements. To pro- 
duce a continuous distribution rather than a finite set 
of levels the length scale associated with this disorder 
must greatly exceed the lattice spacing. Possible mech- 
anisms include a low density of extended defects or a 
density wave which generates a continuum of local envi- 
ronments. The natural culprit would be Bi doping how- 
ever this would appear to be at odds with the fact that 
a single distribution function measured for x = 0.97 ac- 
counts for all specific heat data for x > 0.8. There are 
in fact indications of disorder for both end members of 
the series. While the high T specific heat anomaly for 
Pr 2 Ru 2 07 indicates quasi-long-range Neel order below 
Tjy(Ru) = 165 K, the phase transition is accompanied by 
a difference in the field cooled and zero field cooled sus- 
ceptibility Furthermore, with a residual resistivity of 1.5 
mSlcm the metallic end member, Bi 2 Ru 2 C>7 appears to 
be severely disordered. These results suggest that a com- 
mon low density of large defects principally controlled by 
the ruthenium sub-lattice may influence Pr 2 - x Bi :r Ru2 07 
at least through the metallic part of the phase diagram. 

Pyrochlore systems are indeed known to suffer from 
various types of clandestine disorder. Y 2 Mo 2 07 was for 
example recently shown to have significant disorder on 
the Mo sublattice [J^j. Virtually undetectable through 
diffraction this quenched disorder nonetheless causes spin 
freezing for T < T g = 22.5 K. Just as non-Kramers dou- 
blets are particularly sensitive to disorder, so are highly 



frustrated magnets so that modest levels of disorder that 
might not affect more conventional magnetic materials 
can have qualitative effects in these systems. Materials 
such as Pr 2 _ x Bi a; Ru 2 07 and Pr 2 Ir 2 07 [J3 that arc both 
magnetically frustrated and contain a non-kramers ion 
are particularly susceptible to these effects. 

Non-Kramers magnetic ions include Pr 3+ , Pm + , 
Sm 2+ , Eu 3+ , Tb 3+ , Ho 3+ and Tm 3+ from the rare earth 
series and U 4+ from the actinides. There are a number of 
materials of current interest where distributed lifting of 
a non-Kramers degeneracies might account for unusual 
low— T anomalies. These include nominally stoichiomet- 
ric materials such as Tb 2 Ti 2 7 [Hj, and Ho 2 Ti 2 7 [iH ]. 
There may also be specific interesting effects of dop- 
ing non-kramers systems as in LiH0.jY1_.rF4 [2(| and 
Tb a; Y 2 _ 2 ;Ti 2 07 |21| where symmetry breaking impurities 
may act as a random transverse field. 
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